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ON THE SUPER-EARTHS LOCKED IN THE 3:2 
MEAN-MOTION RESONANCE 

A. Lacny 1 and E. Szuszkiewicz 2 



Abstract. The first study of migration-induced resonances in a pair of 
Earth-like planets has been performed by Papaloizou and Szuszkiewicz 
(|2005[l . They concluded that in the case of disparate masses embedded 
in a disc with the surface density expected for a minimum mass solar 
nebula at 5.2 au, the most likely resonances are ratios of large integers, 
such as 8:7. For equal masses, planets tend to enter into the 2:1 or 3:2 
resonance. In Papaloizou and Szuszkiewicz (|2005p the two low-mass 
planets have masses equal to 4 Earth masses, chosen to mimic the very 
well known example of two pulsar planets which are close to the 3:2 
resonance. That study has stimulated quite a few interesting questions. 
One of them is considered here, namely how the behaviour of the plan- 
ets close to the mean-motion resonance depends on the actual values 
of the masses of the planets. We have chosen a 3:2 commensurability 
and investigated the outcome of an orbital migration in the vicinity of 
this resonance in the case of a pair of equal mass super-Earths, whose 
mass is either 5 or 8 Earth masses. 



1 Introduction and numerical set-up 

In the first study of migration-induced resonances in a pair of low-mass planets 
Papaloizou and Szuszkiewicz (|2005p performed simulations of two interacting plan- 
ets embedded in a disc with which they also interact (see also Szuszkiewicz, this 
volume). The simulations were carried out for a variety of planet masses, initial 
orbital separations and initial surface densities in order to explore the possible 
outcomes. In the case of planets with equal masses when the relative migration 
is slow, Papaloizou and Szuszkiewicz (|2005p have found that the planets become 
trapped in the nearest available resonance. In their study, the two low-mass plan- 
ets each have a mass of 4 Earth masses, chosen to mimic the very well known 
example of two pulsar planets which arc close to the 3:2 resonance. In this work 
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we have modelled the evolution of a pair of planets of equal mass in a gaseous 
protoplanetary disc, again in the vicinity of a 3:2 mean-motion resonance, but for 
two different, larger, values of the mass, namely, 5 and 8 Earth masses. 

We have used the hydrodynamical code NIRVANA, originally written by Zicgler 
(|1997[) and then adapted to disc-planct simulations. More details about the code 
can be found in the work by Nelson et al. (|2000|) . The kinematic viscosity of 
the disk is zero and no magnetic field is present. The accretion of matter onto 
planets is not included. The units used in the code are the following. The mass 
of the central object is taken to be the mass unit. The initial distance r p i of the 
inner planet from the central body is the length unit. The time is measured in 
(GA/ Q /r P 2 3 ) ■ The profile of the initial surface density in the disc is essentially 
flat and it is described in detail in Papaloizou and Szuszkiewicz (|2005[) . The initial 
value of the surface density So at the locations of the planets is Eo = 6.0 x 10~ 4 
in our units. The computational domain in polar coordinates extends from 0.33 
to 4 length units in the radial direction and covers all azimuthal angles from to 
2tt, and was divided into 384 radial slices and 512 angular sectors. The disc model 
is locally isothermal with aspect ratio H/r = 0.05 

2 Results 

In order to investigate how orbital evolution in the vicinity of a 3:2 resonance 
depends on the mass of the planets, we first considered two planets with 5 Earth 
masses each and then we repeated the calculation for two planets with 8 Earth 
masses each. We assume that previous evolution brought the two planets to a con- 
figuration close to the 3:2 commensurability. For the low-mass planets embedded 
in a disc with the surface density distribution given in the form E(r) ~ r~ a one 
can estimate the characteristic time of migration r on a circular orbit from the 
following equation (see Tanaka, Takeuchi and Ward (|2002[) ) 



where M c is the mass of the central body, M p is the mass of the planet, r p is the 
semi-major axis of the orbit of the planet, fl p is the angular velocity of the planet 
and c is the sound speed in the disc at the distance r p from the central body. 
We expect that each pair of planets will enter the 3:2 mean motion resonance 
(MMR) because of their convergent migration, clearly following from Equation 
(|2.1[) . Among two planets with the same mass, the planet located further out 
from the central body will migrate faster. 

To analyse the system dynamics close to the 3:2 commensurability in the pres- 
ence of the disc we monitor the evolution of the orbital elements, mainly the 
semi-major axes, eccentricities and resonant angles defined as follows 




(2.1) 



<f>i = 3Ai — 2A 2 - lu. 



(2.2) 
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Fig. 1. Resonant angles (upper panel), ratio of the semi-major axes (middle panel) and 
orbital eccentricities (lower panel), crosses - the inner planet, squares - the outer planet, 
of the systems of 5M® planets - left column, 8M0 planet - right column. 



where i = 1, 2, Aj is the mean longitude of the i planet and Wi is the longitude of 
the pericentre of the i th planet. To identify the occurrence of the commensurability 
we have used the fact that near the 3:2 resonance, the resonant angles (|2.2[) and 
u>\ — UJ2 should librate about equilibrium values. In both cases the planets of 5M ffi 
and 8M0 entered the 3:2 MMR and stayed there throughout the whole experiment 
from the beginning till the end, for 25000 dimensionless time units as is seen in 
Figure [TJ 
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3 Conclusions 

The planets with 5 and 8M S entered the resonance and continued evolution in this 
state till the end of our experiment. The mean values of the orbital eccentricities 
are excited to the value of ~ 0.008 in both cases and keep oscillating around this 
value. The eccentricities of the internal and external planets are roughly the same. 
The ratio of the orbital semi-major axes slowly approaches the exact resonance. 
A full discussion of these results will appear in a forthcoming paper. 
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